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1. Introduction 

I describe some of our recent results on computing contributions of higher dimension 
operators to effective AS = 1,2 hamiltonians within the Standard Model. These correc- 
tions are due to higher order terms in heavy quark mass expansion {m^ 1 ) and require 
thorough numerical study before being safely neglected. 

2. Corrections to K — ► nn decays: AS = 1 effective hamiltonian. 

The short- distance analysis after removing the iy-boson and the heavy quarks results 
in the effective AS = 1 hamiltonian of the following form |2j 



where Gf is the Fermi constant, V stands for the Cabibbo-Kobayashi-Maskawa mixing 
matrix, r = —V t dV t * s /V U( iV* s , Zi(fj) and yi(fi) are Wilson coefficients, and {Qi\i = 1, ...,6} 
is the basis of local operators of dimension six in mass units (four-quark operators). 
In Eq. ([!]) contributions of electroweak penguin operators || f|] are omitted. For more 
precise comparison of theoretical predictions with experiment the existing analysis cannot 
be considered complete and following points require further investigation: 

• Perturbation theory for kaon nonleptonic decays can be improved by using more 
accurate effective hamiltonian. Within the standard approach only leading terms in 
the inverse masses of heavy particles are kept while a proper account of nonleading 
corrections in the inverse mass of charmed quark || can be important. Further 
corrections to the effective hamiltonian appear because the top quark is heavier 
than other quarks and VT-boson. This results in an incomplete GIM cancellation 
and in appearance of a new operator in the effective hamiltonian || . 

• High precision theoretical estimates for kaon decays stumble at a necessity to cal- 
culate mesonic matrix elements of local four-quark operators entering the effective 
hamiltonian ([!]). The only method of computation entirely based on first princi- 
ples seems to be numerical simulations on the lattice (see, e.g. J7j) though so far 
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even this approach has not given unambiguous and stable results. Several semi- 
phenomenological techniques have been developed and applied for computation of 
matrix elements, e.g. ||, though the precision still need to be essentially improved. 
Recently a regular method to evaluate the mesonic matrix elements has been ex- 
ploited H where the effective hamiltonian is represented in terms of the chiral theory 



variables pi and parameters of the chiral representation are determined via QCD 



sum rules JTT|, [12], [13| for an appropriate three-point Green's function. 



Perturbation theory does not take into account soft light quarks and gluons with 
small virtual momenta. They are entirely hidden in matrix elements of local four- 
quark operators. Well known factorization procedure for evaluation of these matrix 



elements [|14[] accounts only for the " factorizable" part of the interaction [p~5] , [16], [T7|j . 
" Unfactorizable" contributions, for example, those corresponding to annihilation 
of a quark pair from the four-quark operator into soft gluons [|T^] are omitted. 
The calculation of these contributions and the generalization of the matrix element 
estimates beyond the factorization framework can be systematically done within the 
approach of ref. JjlJ. 

After decoupling heavy particles (VT-boson, t-, b-, and c-quarks) from the light sector of 
the theory, Eq. (|l|) corresponds to the leading order in inverse masses of these particles. 
The removal of the c-quark however is not very reliable and, in general, requires a special 
investigation. It is not heavy enough in comparison with a characteristic mass scale in the 
sector of light u-, d-, and s-quarks, for example, with the p-meson mass. The nonleading 
terms in the l/m c expansion can, therefore, be important and require a quantitative 
consideration. To compute corrections of order l/m c , the tree level hamiltonian before 
decoupling of the c-quark is used. It reads 

H£ s =i = ^V ud V: s {Qt - (1 - r)Ql) + h.c. (2) 
where Q\ = ^s^^^X^^g^), Ql(r) stands for left(right) handed quark. Performing the 
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OPE and restricting oneself to the first order terms in a s and m c one finds 

H AS=1 = H^ + H^. (3) 

The first addendum in Eq. (|3[) corresponds to leading contributions in l/m c and 
coincides with Eq. ([!]). Second addendum in Eq. (|3|) is the l/m c correction || 

# (8) = - Q? ) + £ cPm s QP) + ^ (4) 

where a basis {Qj 8 ^ = 1, 7} |i = 1, ...,4}) of the local operators with dimension 

eight (seven) in mass units is chosen in the form 

Q2 = ig s SL{JfilaG a ^ - JaG^J^dl, 

Q3 — S l{PaG paflvG V fi + 1uG ufl G ^a-PO^L) 

(8) " a 

= "iSL{G flu a fll j'-f a G a [3Pf3 — PplaGapG ^pU ^ v )dL, 

Qf = s L (D 2 J)d L , = is L {DG vlx G vlx - G u ,DG u ,)d L , 

Q^p — sr{G pvG a p)dL, Q2 = sr(G fj, u G ut j)dL, 

Q? = is R (G ua G afl a Uf ,)d L , Q { p = s R {J^ + P^J^di. (5) 

Here P M = id^ + g s A^ is the momentum operator in the presence of external field = 
Ap: a , t a are generators of the color group SU(3), G^ v = G^ u t a is the gluon field strength 
tensor, J ^ = J2 q =u,d,s(^^ a( l) ta ^ and °> = ^7/*, 7*]/ 2 - In derivation of Eq. (|) u- and 
ci-quark are considered massless, and the first order in strange quark mass was kept with 
the use of equations of motion sP = m s s, Pd = 0, [P M , G^ v ] = iD^G^ u = —ig s J u . 
Straightforward calculation gives the following values for coefficients Q to the leading 
order in a s 

Cf = -2Cf ] = 4Cf = -8Cf } = 2Cf ] = 8Cf = -g, Cf = 0, 



CP = CP = -1 CP = 6 -, cf = 0. (6) 

This generalizes the effective hamiltonian for AS = 1 decays beyond the leading order in 
l/m c expansion. 

To complete our treatment of the local effective hamiltonian we consider the case of 
a heavy top quark. GIM cancellation is not perfect and the quark-gluon operator m s Q^ 5 ' 
appears in the effective hamiltonian already in the first order in a s . This additional 
contribution || reads 

A#(6) = ^VrtVZTC®fr)m.Q®(jM), 

1 Tfi^ 
C {5 \ri = 77 r^(F(x c )-F(x t ))r ] (ri, x q - 1 



11/5 39 \ 

F(x q ) = - _ i^-x, 4 - 7x q 3 + -x q 2 - 19x q + 4 - 9x q 2 \nx q j 



F(x c ) ~ F(0) = | (7) 



The renormalization group factor has the form 



_ , x \ 7 (5) /2& /- / A \7 (5) /2/3 4 / — / \ \7 (5) /2/9 3 



\a s (M w ) J \a s (m b ) J \a s (m c ) J 

where 7^) = —28/3 is the anomalous dimension of the operator m s Q^ [|2U[| , /3 nf — 
11 — |n/, n/ is the number of active quarks flavors. 

Contributions of u- and c-quarks to the real part of the effective hamiltonian cancel 
each other via GIM mechanism and the operator m s Q^ 5 ' contributes to imaginary parts 
of amplitudes and, therefore, can be important in the analysis of direct CP violation. 
However, its contribution is suppressed numerically because rj(fi) < 1 and the function 
F(x) changes slowly. Indeed, at the point Aqcd = 0.3 GeV, fj, — 1 GeV, m t = 130 GeV, 
one has = 0.001, while the numerical value of the Wilson coefficient of the dominant 
penguin operator Q 6 = -8E, =uAs (sL9«)(9ii4) is y& = 0.1. 

Thus, the complete form of the effective hamiltonian up to the first order in m s , l/m c 
and a s with m t ~ M\y is now available. 
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As an example of using the above hamiltonian we consider K — > tttt decays. For this 
end we have to extract an information about the matrix elements of the local operators 
between mesonic states. We start with operators m s Q^ and m s Q^ . These operators 
contain explicitly the strange quark mass and, therefore, in the leading order of the chiral 



expansion they correspond to tadpole terms in the chiral weak lagrangian |2T| that do 
not generate any observable effect and can be neglected in the leading order of chiral 
symmetry breaking p2, p3fl . 



Thus the problem is reduced to the estimation of the matrix elements of the oper- 
ators Ql . At present there is no regular method to calculate them within QCD ex- 
cept direct simulations on the lattice. To estimate at least the scale of nonleading l/m c 
corrections we work with a simplified model that uses factorization. One selects opera- 
tors containing scalar quark currents which can be written as (slG ^cr Atl/ g fi )(g J j(i L ) and 
[s lQr) (qrG ^ u a ^ydi) ■ This step seems to be justified because in the case of dimension six 
operators the similar "penguin- like" structures are strongly enhanced and dominate the 
others. The last simplification consists in the substitution qg s G ^q — > m^qq where ml 
determines the scale of nonlocality of the quark condensate and is defined by the equation 
(qgsG^a^q) = m 2 (qq), mg(l GeV) = 0.8±0.2 GeV 2 @, |2jJ. This substitution is valid 



in the chiral limit for the operator qg s G ^a^q. We suppose that it is justified also in our 
case at least for estimates of the order of magnitude. All above assumptions about the 
factorization procedure in the case of dimension six operators become exact within the 



many color limit of QCD, N c — > oo, to the leading order in N c [26 



(8) 

Thus, only operator Q\ has a nonvanishing matrix element that reads 

(nn\QP\K) = ^(nn\Q 6 \K). (8) 

(8) 

We should note that the coefficients CI are finite to the leading order in a s and indepen- 
dent of renormalization scheme. We can therefore use the leading order values of mesonic 
matrix elements that is consistent up to the considered level of accuracy. The next-to- 
leading a s corrections to Wilson coefficients depend on the renormalization scheme and 
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matching between them and mesonic matrix elements is necessary to make physical am- 
plitudes independent of the renormalization scheme. Thus, account for the first order 
1 / m c corrections results in the shift of coefficients in front of the penguin operator Qq 

7 - (y + 013 m ° r { A 7, - (n, ^ ^ r (8 A (Q) 

+ j, y8 - > ^ 6 ___c 4 j. (9) 

Using the numerical values z 6 = —0.015, y 6 = —0.102 at the point A QCD = 0.3 GeV, 
\i = 1 GeV, nit — 130 GeV one finds numerically relative corrections to the Wilson 
coefficients to be 

z 6 -> ze(l - 0.1), y 6 ^ y 6 (l + 0.01). 

The main correction appears in the real part of the Wilson coefficient of the penguin 
operator Qq. Parametrically, the contribution of dimension eight operators can be as large 
as 50% {rnl/m 2 c ~ 0.5) of the leading term. However much smaller value has been found 
within the simplest factorization framework for meson matrix elements of nonleading 
1 / m c contributions to the kaon decay amplitudes. Large violation of the factorization for 
matrix elements of dimension eight operators seems to be likely and the numerical value 
of nonleading l/m c correction can be estimated only when a self-consistent method to 
calculate these matrix elements within QCD will be available. 

3. Corrections to K° — K° mixing: AS = 2 effective lagrangian. 

The effective local AS = 2 lagrangian for the K° — K° mixing in the leading approx- 



imation in the heavy charmed quark mass m c is well known p7j , p8| . The corresponding 
hadronic matrix element of the effective local Las=2 lagrangian between K° and K° states 
(K c \k')\L &s=2\K z '(&)) has been intensively studied during several last years with different 



techniques {e.g. W7L IT9], RO, RTl 132]). However it has been pointed out in 133 that 



the local effective hamiltonian does not exhaust the physics of AS = 2 transitions. It 
cannot account for the long distance contribution which is present in the initial Green's 
function for the matrix element of the K° — K° mixing and is connected with the propa- 
gation of the light w-quark round the loop of the box diagram. This contribution is purely 
nonperturbative and ultimately depends on infrared properties of QCD. 



Nevertheless there is one point which can be essentially improved just within pertur- 
bation theory for the standard model. It consists in the calculation of corrections in the 



inverse mass of charmed quark to the local part of the effective lagrangian p4 |. These 
corrections are represented by local operators with dimension eight in mass units. 
Because the AS = 1 lagrangian has the form 

Las=i = ^J,J^ ( 10 ) 

where J M = Ql^uVql is the weak charged hadronic current, Q = (u,c,t) T , q = (d,s,b) T , 
the matrix element M of the transition is represented by 

out (K°(k')\K°(k)) in = i(2ix) A 5{k - k')M, 

M= l -J dx(K°(k')\TL AS=1 (x)L AS=1 (0)\K (k)). (11) 



Eqs. (|T0| - |TT|) are valid for the t-quark much lighter than the jy-boson. This is not the case 
anymore but for our purpose it is inessential and in the following we neglect the t-quark 
admixture and restrict ourselves to the simplified model with two generations. 
The effective AS = 2 lagrangian can be written in the form 

cos 9 C \ 2 

v2 

where 9 C is the Cabibbo angle. Here 



/ AGp sin 9 C cos 9 C . 
L AS =2= — ^ -) (L H + L L ) (12) 



L H = i J T H (x)dx, T H = T cc - T cu - T u 



is the heavy part of the whole effective AS = 2 lagrangian containing loops with virtual 
heavy c-quark in the intermediate state, while 

L L = i J T L (x)dx, T L = T uu 

describes the light part of the transition. We introduced useful notations 

T cu (x) = Ts L j a u L c L 'j a d L (x)s L 'yi3C L u L 'ji 3 d L (0) 
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and so on. Both Lh and Ll separately require some regularization because they are 
ultraviolet divergent. The dimensional regularization is not convenient in this case due 
to the presence of the 75-matrix. The Pauli-Willars regularization introduces a regulator 
mass that makes difficult to perform explicit calculations. We will use the regularization 
which is free of these shortcomings in our particular case. Namely, let me define the 
regularized quantities L HL by the equation 



L H ,L = * J T H!L (x)(-fx 2 x 2 ) e dx 



where e is a regularization parameter and \x represents the mass scale analogous to one of 
dimensional regularization. 

Now for the heavy part of the effective AS = 2 lagrangian Lh we develop a regular 
expansion in the inverse charmed quark mass in the following form (from now on we omit 
the index "R") 

16n 2 L H = C Q (ii, m c )O (p) + ^ Cj(fi, m c )Oj{ij) 

j 

where Cj are coefficient functions depending on the heavy quark mass m c and Oj are 
the local operators built from the light (u, d, s) quark fields only. If we split the whole 
lagrangian into the sum 

L H = LP + Lg } (13) 
then 167r 2 L^ = —m 2 {sL^adL) 2 is the well known result of Gaillard and Lee |14]]. The 



rest part of Eq. (|T^) contains local operators which have dimension eight in mass units. 
A convenient form of the operator basis is 

Op = SLladLSLliiF^dp, A = s L ~i^D u) d L s L -i {fl D v) d L , 

B = SLj^D^dLSL^D^dL, 

Oc = SLjadLSL^D^Do, + -D Q 7 M -D M )rf L s —^ — —{s L -i a d L ) 2 

where J^D^ = {pf^D u + •y u D^)/2. The direct calculation gives 



where C = 0.577... is the Euler constant. After performing a renormalization (say, a 
minimal subtraction of the pole term) we will have the finite quantity and the parameter 
\i recalls the necessity to have the proper short distance contribution of w-quark. In this 
order of the expansion in m" 1 the dependence on this parameter is explicit contrary to 
the leading order that is finite and does not depend on \x. The heavy and light parts of 
the whole lagrangian must be defined simultaneously in the coordinated way. The pole 



part of Eq. (|lj) is cancelled by the corresponding divergences of the light part due to 
GIM mechanism. Consider now the light part. The operator product expansion for the 
amplitude T L [x) in x 2 at x 2 — > has the form 

T L (x) = -^{s Lla d L f + (20 p + 20 A + B + Oc) • (15) 

The short distance contribution of the light part does cancel divergences of the heavy 
part. More technically we extract the short distance contribution of the light part by 
splitting the entire light part into a sum 



L L =i J T L (x){-fi 2 x 2 fdx = i JT L {x){f{x,x) + f{x,x)){-fi 2 x 2 ) e dx = L s L H + L 1 [ G (16) 
where 

L S L H = ij T L (x)f(x,x)(-fi 2 x 2 ) e dx, L L L G = ij T L (x)f(x,x)(-fi 2 x 2 ) e dx, 

and f(x, x) + f(x, x) = 1. The smooth generalization of the step ^-function is chosen for 
functions f(x,x) and f(x,x) 



x 2n Tl _ N (-x 2 ) n 



f(x, X) — _ow> f{ X i X ) 



x 2n + (-x 2 ) n x 2n + (-x 2 ) n 

in such a way that the function f(x,x) cuts out the short distances only (up to x) and 
the function f(x,x) does the long distances. The short distance contribution of the light 
part is 

16vr 2 Lf = {s Lla d L f n,n 4 + \ ( 2 °f + 20a + O b + O c ) (- + ln/i 2 x 2 ). 

sm(7r/n)x^ 3 e 
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Now the whole answer is 

IQti\L h + L L ) = 1Gtc 2 L l l g + {s Lla d L f ( -m\ + " / " ' 



sin(7r/n) x" 1 



2 / 4e~ 2C 4\ 2 

- - (20 p + 20 A + B + Oc) f M^pf) + 3 J - (17) 

The long distance part of the whole light lagrangian cannot be calculated due to 
strong infrared problems and requires some low-energy model, lattice or chiral effective 
theory |35| . Numerical estimates for corrections depend on kaon-antikaon matrix elements 



of the local operators Op, — Oc for which factorization was used 

(K°(k)\(s L -f a d L ) 2 \K°(k)) fact = (1 + ±-) = Z, 

(K°(k)\Op\K°(k)y act = -5 2 f^y^) , (K (k')\O A \K°(k)) fact = -5 2 -^ (^^j 

where the parameter 5 2 is defined by the relation |36| {Q\sL^^F^ a dL\K°{k)) = —ik^f K 5 2 
and all other matrix elements vanish. From Eq. ( |TTD one finds 

WA L H + Li-) = Z(- m l + * + ^(ln(g) + \) + J*"). 

Let me estimate x 2 . Eq. ([15|) gives 

r *(*) = 4^(l-^ + (**)). (18) 

Then at <5 2 x 2 = 3 where the expansion ([18]) blows up and at n = 00 we obtain the following 
representation for nonleading corrections 

2r7 / 4 5 2 4 5 2 ,4<5 2 , n ^ 35 A 
-m\Z 1 - -— - -—In — -2C + -)). 
\ dm* 3 m A c dm* 24 y 

Numerically, at 5 2 = mg/4 we get 

-m\Z{\ - 0.2 + 0.3) = -m\Z{\ + 0.1). 

Defining the value of x 2 by the relation 5 2 x 2 = 1 we get 

-m\Z{\ - 0.5 + 0.1) = -m 2 Z(l - 0.4). 
11 



4. Conclusions. 

The corrections in the inverse mass of the charmed quark to the effective AS = 1,2 
lagrangians can be important numerically for precious comparison theory with experi- 
ment. The AS = 2 lagrangian reveals at this level an explicit dependence on the boarder 
between short and long distances. The OPE at x 2 — > for the light part of the lagrangian 
is analyzed for determining the convergence scale. Within vacuum dominance approxi- 
mation this scale is large enough and is given by the parameter 5. The corrections can 
change the transition matrix element considerably depending on the concrete choice of 
the scale x. This means that the long distance hadronic contribution is important to 
stabilize the result in the region where the expansion ( |i"8|) fails, in other words the inte- 
gral I/£ changes quickly for x 2 in the region 1/5 2 < x 2 < 3/5 2 in order to compensate 
corresponding changes of the short distance part (|T5|). 
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